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Abstract 
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1.  Introduction 


Assume  y  to  be  a  bounded  measure  on  Let  f  be  a  kernel  changing  y 

into  the  transformed  function 

CO 

(1)  g(t)  :=  /  f  (t  ,u)djJ(u) . 

0 

In  the  literature  dealing  with  operational  calculus  the  recovery  of  y 
from  g  is  mostly  accomplished  by  complex  inversion  formulas.  One  knows  that 
for  a  variety  of  integral  transforms  (i.e.  Laplace,  Stieltjes),  real  inver¬ 
sion  formulas  are  available  [6,10,12].  These  formulas  can  also  be  obtained 
as  applications  of  classical  results  from  probability  theory  such  as  the 
weak  law  of  large  numbers  and  the  central  limit  theorem. 

To  abbreviate  the  writing  we  put  for  0<y1<y2<o° 

(2)  v(yl>y2}  =  ^{yi}  +  ii^i*y2)  +  ^{y2} 


•Access' ;  r.  V 

’MIS 

V>TIC  1  -h 
Jnunro.!t»o-> 5 

iustirica:. . 

•  l3tributt 
a  va '.lab  ill  tv 

[Avail  - 


2 .  Inversions  using  lattice  variables 

Let  X  be  a  non-negative  lattice  variable  with  mean  u  >  0  and  finite 
2 

variance  o  >0.  Take  X^,X2>...  a  sample  from  X.  Define  the  sum 

S  (“)  *  X.  +  X_  +...+  X  .  n  >  1 
n  1  2  n 

and  put 

d^(u,t)  =  PtS^fu)  s  [nt J } . 

p 

By  the  weak  law  of  large  numbers  we  know  that  Sn(u)/n  -*■  u  while  [nt]/n  ■*  t. 
Hence  dn(u,t)  +  0  if  t  <  u,  while  d^Cu.t)  -*■  1  if  t  >  u.  For  u  =  t  note  that 
{ [nu]  -  nu}//n  0  while 


S  (u)  -  nu 
(u>u)  =  P{_Q - 

a/n 


[nu] 


nu 


aSn 


}  -4 


by  the  central  limit  theorem  and  the  uniform  convergence  to  the  normal  law 
[7,  p.  139].  Hence 


dn(u,t) 


0  t  <  u 


1_ 

2 


t  =  u 


1  l  t  >  u  . 


By  Lebesgue's  theorem  and  (2)  for  0  <  y^  <  y^  <  00 

OO 

(3)  yty.jy.}  =  lim  /{dn(u,y2)  -  dn(u,y1) }dp(u) 

n-*»  0 


is  a  potential  inversion  formula. 

Moreover  by  a  local  limit  theorem  (4,  p.  225]  or  [1,  p.  233] 


(4) 

I 


hn(u,t)  :=  o/Zrrn  PiS^fu)  =  [nt]} 


0 

1 


t  ^  u 
t  =  u  . 


By  proper  approximation  (4)  will  suggest  a  function  of  the  form 


e  (v)  =  ch  (v,l) 
n  n 


0 

1 


v  /  1 
v  =  1 


for  a  constant  c;  this  implies  an  inversion  formula  for  the  point  mass 
y,  i.e. 


(5)  p{y}  =  lira  /  e  Ady(y). 

n-*»  0  y 

Before  turning  to  examples  let  us  introduce  the  abbreviation 
(0  <  yl  <  y2  < 

(6)  An  =  {m:  [ny^  +  1  <  m  <  [ny2]>. 

Example  2.1.  Poisson  distribution  and  Laplace  transform. 

Put  for  the  Laplace  transform  of  p, 

OO 

(7)  g(t)  :=  /  e'tudp(u). 

0 

Then  for  m  e 

(8)  (-l)ng(n)(t)  =  /  e'tUundp(u). 

0 


Let  X  be  Poisson  with  mean  u  so  that 


d 

n 


(u,t) 


li]* 

ra=0 


-nu 


(nu)m 

m! 


By  (2)  and  (3)  we  get 


oo  m 

P{y1iy2^  =  lira  I  /  e'nu  y  ■  dy(u) 


n-+°°  meA  0 
n 


=  lim  £ 
n-*°°  meA 


mi 


n 


This  inversion  formula  can  be  found  in  [12,  p.  295].  It  can  be  fruitfully 
used  in  giving  a  probabilistic  proof  of  Bernstein's  theorem  on  completely 
monotone  functions;  see  [2,  p.  191]. 

For  the  point  mass  we  get  from  (4) 

,  ,  ^  k -  -nu  (nu) 

hn(u,t)  =  e  -frt]!  -  • 


Stirling's  formula  easily  yields  c  =  1  and 
efi(v)  =  exp  -  n{v  -  1  ♦  log  v). 

Collecting  results  we  get 

Theorem  2.1.  Let  p  be  a  bounded  measure  and  g(t)  its  Laplace  transform. 
Then  for  0  <  y^  <  y ^  <  00 


y{y 


ry2} 


lira  l 
n-*30  meA 


(-nu) 

m! 


m 


n 


and  for  0  <  y  < 


OO 


2. 


p{y}  =  lim(-  ^)ng(n)(^). 


n-wo 


The  latter  formula  occurs  also  in  [12,  p.  298] 


Example  2.2.  Pascal  distribution  and  Stieltjes  transform. 

In  the  sequel  we  also  deal  with  the  ^generalized)  Stieltjes  transform. 

OO 

(9)  4>(t)  :=  /  (u  +  t)"Pdp(u) 

F  0 

where  p  >  0. 

In  evaluating  successive  derivatives  of  the  next  lemma  will  be 
useful . 

Lemma  2.2.  For  a,  ecu,  -p  ^  JJq,  u>0,  t>0 


(10)  DaltbDC(u  .  t)-P)  .  C)(-Da,Ctb(u 


+  tl 


-p-a-c 


F ( -a . -h . -t 


in  particular 

(11)  Da{tP*a*C-1DC(u  ♦  t)-P)  =  f.(f  *r»  *  C)(-l)CtP<C-1ua(u  *  t)-P-3-C. 

Proof.  Apply  Leibniz's  rule  for  differentiation  of  products;  and  then,  for 

a  en 

a  (a)  (6) 

F(-a,-8,-y;x)  -  l  -^-(-x)ra  . 

m=°  Cm) 

To  prove  (11),  note  that  F(-a,y,Y;x)  =  (1  -  x)a  if  a  eW  by  [5,  p.  1040] .  Q 


Let  X  be  geometric  with  mean  u.  Then 


2.5 


[nt]  . 

,  ,  . ,  V  .n+m-1.  m,,  .  n 

d  (u,t)  =  l  (  )P  (1  -  p) 


where  u  =  EX  =  or  P  =  (u  +  1)  *.  For  the  point  mass  we  note  that 


1  m 

h  (u,t)  =  o/2irn(n+m  )  —  u-~ — - 
n  m  (l+u)n+m 


where  m  =  [nt] .  Stirling's  formula  again  applies  to  yield 


,  ,  ,2i/v.2n 

n  Vl+v' 


Theorem  2.5.  Let  y  be  a  bounded  measure  and  <p ^ ( t )  its  (ordinary)  Stieltjes 
transform.  Then  for  0  <  y^  <  y^  <  00 


■ li,n  nsy  1 

n->oo  me  A 


for  o  <  y  <  °°, 


y{y)  =  lim  -?^yf)n{y2n'1(-D)n"1<|)1(y)}. 


Proof.  From  (3)  and  the  current  choice  of  d  we  find 
— -  n 


-  u.  i  ("T1)/ 


,  00  n 

n+m-K  r  u 


dp(u) 


n-*°  me  A  0  (u+1) 


Put  a=n,  b=n+m-l,c=m-l,p=lin  (11)  to  find  the  first 


2.6 


formula.  To  get  the  second,  put  a  =  n  +  m  -  1,  b=n,  c=0,  p=l  in  (10) 
and  -ecall  that  F(a,b,c;x)  =  F(b,a,c;x). 

For  the  point  mass  we  get  from  (5)  and  (12) 


oo  n 

y{y}  =  lim(4y)n/  — - — =—  dy(u). 
n-»»  0  (u+y) 

So  we  take  a=n,b=2n-l,  c  =  n  -  1  and  p  =  1  in  (11). 
Example  2.3.  Binomial  distribution  and  Stieltjes  transform. 
Let  X  be  Bernoulli  on  0  and  1  so  that 


P{X  =  k)  = 


k  =  0 
k  =  1. 


Then 


P(Sn(v) 


[nt]} 


-  [f 1 

1=0 


r  n.  ^ . 
(^)v  (1  - 


v) 


n  -Z 


suggests  a  potential  inversion.  However  we  should  adapt  it  so  that  v  e F 
rather  than  v  e  [0.1],  So  define 


d 

n 


(u,t) 


l  (-)u  (1  ♦  U)  . 

1=0  C 


Theorem  2.4.  Let  p  be  a  bounded  measure  and  ^(t)  its  (ordinary)  Stieltjes 
transform.  Then  for  0  <  y^  <  y^  <  00 


where 


y{y1;y2l  =  ii»  I  (tj) 

n-*“  £eB  ^  J 


=  lim  l  THrnT(J)Dn'1{tn_£4.  (t)}  , 

n-+°°  leB  {n  1J  *  1  1  1  1 

m 


Bm  =  r  1J  <  £  < 


Ly2+i 


1 


The  proof  is  basically  the  same  as  that  of  the  previous  theorem. 


Inversions  using  also  densities 


3 . 

Let  be  a  sequence  of  independent  non-negative  random  variables 

with  common  distribution  F,  mean  1  and  finite  variance.  Let  also  Y  ,Y7>... 
be  such  a  sequence  but  with  distribution  G,  and  independent  of  the  sequence 

{X.  }. 

1 

Define  S  =  X,  +  ...+  X  and  T  =  Y  +...+  Y  and  look  at 
n  1  n  n  1  n 

d  (u)  =  P{S  /T  <  u}; 
n  n  n 


by  combining  a  weak  law  and  the  central  limit  theorem,  one  easily  finds  that 

rO  u  <  1 


00 

d  (u)  =  /  F(n)(ut)dG(n)(t)  - 
* 1 


u  =  1 
u  >  1 


where  F  ^  (G  )  is  the  n-fold  convolution  of  F(G)  with  itself.  The  particular 
case  where  Y^  =  1  is  important  and  yields  d^(u)  =  F^fnu). 

The  potential  inversion  formula  reads  now 

00 

(13)  u{y1;y2>  =  lim  /{dn(— )  -  dn(^-)}dy(u) 

n-*°°  0  yl  y2 

while  a  point  mass  formula  can  be  derived  by  using 

e  (v)  =  —  vd1 (v) 

n  /R  n 

if  d  is  absolutely  continuous, 
n 

Lxamplc  3.1.  Lxponcntial  distribution  and  Laplace  transform. 

Here  Y.  =  1  and  X^  has  an  exponential  density  with  mean  1.  Hence 


and  so  we  have 


Theorem  3.1.  Let  p  be  a  bounded  measure  and  g(t)  its  Laplace  transform. 


for  0  <  y  <  y_  < 


n  gl  '  Cjr) 


,  ,  )2  (-n)n  s  't 

"W  •  22  I  -Tw—1 


Proof.  From  (13)  and  (14)  we  see  that 


y{y.;y  }  =  lim  /dy(u)  f  *  e"vvn_1dv 

n-><»  0  nuy2  1  ’ 


oo  n 


lim  '  "TtnT  ~nTT  /  e‘u^(-u)ndp(u) 
n-*»  yj  1  1  t  0 


by  Fubini’s  theorem. 


The  point  mass  formula  coincides  with  that  of  th.  2.1. 
The  appealing  form  of  (15)  suggests  that  if  the  sequence 


(n)  ,n. 
(_n,n  J(y) 

k  (t)  :=  ~ 

n  f(n)  tn+l 


converges  a.e.  on  ]R  +  and  is  bounded  by  an  L^- function  then  y  is  absolutely 
continuous  and  its  derivative  is  given  by  the  inversion 


3.3 


The  inversion  formula  (16)  is  well  known  [6,12];  its  probabilistic  proof  due 
to  Feller  [3]  has  been  the  inspiring  source  for  the  current  paper.  Formula  (16) 
has  been  used  by  Jagerman  [9]  to  perform  numerical  inversion  of  the  Laplace 
transform.  For  other  applications  see  a  paper  by  Vinogradov  [13]. 

Example  3.2.  Exponential  distributions  and  Stieltjes  transform. 

Let  F(x)  =  G(x)  =  1  -  e“X;  then  dn(u)  turns  out  to  be  a  beta  distribution 


on  K  ,  i.e. 
+ 


P  ,  u  n-1 

d  (")  .  J  -i - ^dv. 

n  f2(n)  0  (1+V)zn 


Theorem  3.2.  Let  p  be  a  bounded  neasure  and  <]>.(t)  be  its  (ordinary)  Stieltjes 
transform.  Then  for  0  <  y1  <  <  00 


u{y,;y9}  =  lim  /2  Dn{t2n'1(-D)n'1<J>1  (t)  }dt 

1  2  n^°  r( n)  y. 


n-x»  f^(n)  y ^ 


for  0  <  y  <  00 


p{y}  =  lim  D20  1  Cy)  }  ■ 


Proof.  We  easily  get 


u{y1;y2} 


lim  T-^-  /2  tn_1dt /  U  ■  2  -  dp(u) 
n-x»  r(n)  y2  0  (u+t) 


*  .  •  ,*•  -  V  *  »'*  •> '*  •**  < 


We  get  the  two  formulas  by  applying  lemma  2.2  with  the  two  choices 


a=n,b=2n-l,  c=n-l,  p=l  in  (11) 
a  =  2n  -  1,  b  =  n,  c  =  0,  p  =  1  in  (10) 
as  in  the  proof  of  th.  2.3. 

For  the  point  mass  we  get  en  as  in  (12);  the  resulting  formula  comes  from 
the  second  choice  above.  □ 

For  the  above  formulas  see  [12,  p.  350].  For  applications,  see  [8). 

The  last  example  shows  how  an  almost  trivial  probabilistic  argument  widens 
the  applicability  of  the  previous  examples. 

Example  3.3.  Exponential  distributions  and  generalized  Stieltjes  transforms. 

Return  to  the  previous  example  but  look  at  S  ,  and  T  instead  of  S 

r  n+b  n+a+p  n 

and  where  a,b  e  2Z  are  fixed.  Then 

S  ,  S 

lim  P{^r~ ~ —  <  u}  =  lim  P{^-  <  u}. 
n^°°  n+a+p  n-+°°  n 


Hence  we  have  another  candidate  for  d  ,  i.e. 

n 


d 

n 


(u) 


T  (2n+a+b+p) 
r(n+a+p)r  (n+b) 


u  n+b-1 

f  _ X _ 

,,  .2n+a+b+p 

0  (1+v)  r 


By  the  usual  procedure  we  obtain  also 

b-a-p 

,  .  2  r2/v  ,2n+a+b+p 

=  v  1.7777}  v 


If  we  now  go  through  the  same  routine  as  in  example  3.2  with  a  replaced 
by  n  +  b  and  c  by  n  +  a  we  obtain  the  following  results  of  E.R.  Love  and 
A.  Byrne  [10,11].  We  note  that  we  have  to  restrict  p  to  N  in  view  of  the 


3.5 


probabilistic  set  up. 

Theorem  3.3.  Let  u  be  a  bounded  measure  ;  let  be  its  generalized  Stieltjes 

transform.  Let  a,b  e  7L  .  Then  for  0  <  y^  <  <  00 


VV  ■  £  rJU  {  °”*b{b“P -‘(-d, -♦„»,}  d„ 


for  0  <  y  <  00 


r  i  22n+a+b+pr(p)  r>n+br  2n+a+b+p-l,  „,n+aA  ,  ^ 

M{y}  =  lim  -Tr-  — yp  {y  C-D)  *(y)}. 


n-x» 


T(2n  +a+b+p) 


Note  that  by  using  different  random  variables  one  can  get  a  variety  of  real 
inversion  formulas  for  the  same  Stieltjes  transform  <j>j .  Needless  to  point 
out  how  th.  3.2  follows  from  the  last  result. 

I 

I 

i 

I 


i 


i 


i 


4.  Remarks. 


1.  The  main  drawback  of  the  method  used  in  this  paper  is  that  the  relevant 
integral  transform  has  to  be  recovered  from  a  potential  inversion  formula. 

2.  The  main  virtue  of  the  paper  lies  in  singling  out  the  relevance  of  the 
asymptotic  behavior  of  the  sequence  d^  for  n  -*■  probabilistic  arguments  often 
provide  this  behavior  immediately .v 

3.  There  exist  additive  counterparts  for  the  formula  (13).  For  if 


d 

n 


(v) 


0 

1 

'  2 

1 


v  <  0 
v  =  0 
v  >  0 


where  v  e  p,  then  a  potential  inversion  formula  for  a  transform  on  F  is 


y(y,  ;y9>  =  n™  /{d  (u  -  y,)  -  d  (u  -  yj  }du(u) . 
1  2  n-°R  n  1  n  2 


1  11  si  u 

The  best  known  example  is  d  (v)  =  —  /  -  du  leading  to  the  inversion 

_oo 

formula  for  the  Fourier  transform 

Y  n 

M^l ;y2>  =  lim  /2dx{|^-  /  e'lzx<})(z)dz} 
n-*“  y^  -n 


where  4>(z)  =  /e  dp(u) 


v^Tv  4x2 


— *  ^  >  in  ^ 

Another  example  is  d  (v)  =  -  /  e  dx  leading  to  a  real 

n  -co 

formula  for  the  Gauss-Weierstrasz  transform  [6] 


inversion 


»x(t) 


_  CO 

*2?  ^ 


exp  -  |(u  -  t.)2dy(u). 
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